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Quantum metrology overcomes standard precision limits by exploiting collective quantum super¬ 
positions of physical systems used for sensing, with the prominent example of nonclassical multi¬ 
photon states improving interferometric techniques. Practical quantum-enhanced interferometry is 
however vulnerable to imperfections such as partial distinguishability of interfering photons. Here 
we introduce a method where appropriate design of the modal structure of input photons can allevi¬ 
ate deleterious effects caused by another, experimentally inaccessible degree of freedom. This result 
is accompanied by a laboratory demonstration that a suitable choice of spatial modes combined 
with position-resolved coincidence detection restores entanglement-enhanced precision in the full 
operating range of a realistic two-photon Mach-Zehnder interferometer, remarkably around a point 
which otherwise does not even attain the shot-noise limit due to presence of residual distinguishing 
information in the spectral degree of freedom. Our method highlights the potential of engineering 
multimode physical systems in metrologic applications. 


Quantum phenomena can facilitate and boost the per¬ 
formance of imaging techniques HHa , sensitive measure¬ 
ments in delicate materials 0 0, as well as detection 
schemes probing subtle physical effects such as gravita¬ 
tional waves [7]- These strategies rely on preparing col¬ 
lective superposition states of multiple probes (for ex¬ 
ample photons, atoms) to achieve precision enhancement 
beyond standard limits IHHII]. In the optical domain, 
a common strategy for collective state preparation is to 
realize multiphoton interference in linear circuits, for ex¬ 
ample free space or integrated interferometers [MU], 
fed with nonclassical states of light. Attainable preci¬ 
sion can be however dramatically vulnerable to residual 
distinguishing information between interfering photons. 
Standard methods to improve indistinguishability based 
on filtering are often inadequate, in particular introduc¬ 
ing attenuation which may easily diminish the overall 
benefit of collective state preparation. 

The purpose of this paper is to analyze the interplay 
between degrees of freedom with different experimen¬ 
tal accessibility in two-photon interferometry, which is a 
canonical example of a quantum-enhanced measurement 
[22]. We demonstrate that detrimental effects caused by 
distinguishing information present in one degree of free¬ 
dom that is beyond experimental control or lacks techni¬ 
cal means to improve indistinguishability, can be allevi¬ 
ated by mode engineering in another degree of freedom, 
even though these two remain completely uncorrelated. 
This feature is investigated in the case of local phase es¬ 
timation, whose precision becomes strongly dependent 
on the operating point if the two photons feeding the 
interferometer exhibit residual distinguishability. It is 
shown that a carefully designed preparation and detec¬ 
tion scheme for a degree of freedom other than the one 
causing distinguishability allows one to restore quantum- 
enhanced precision in the entire operating range of the in¬ 
terferometer. We attribute this effect to non-trivial com¬ 
bination of one- and two-photon interference that turns 


out to augment phase sensitivity beyond the shot-noise 
limit. We also present an example indicating that simi¬ 
lar enhancement occurs also at higher photon numbers. 
Prospectively, the results reported here may provide an¬ 
other class of strategies to mitigate effects of imper¬ 
fections and environmental noise in quantum-enhanced 
metrology [23H29] . 

The theoretical analysis is complemented with an ex¬ 
periment investigating sensitivity of a balanced Mach- 
Zehnder interferometer fed with photon pairs. We de¬ 
termine the precision of local phase estimation around 
the operating point when the photons coalesce pairwise 
at the interferometer output ports. In this regime, the 
residual spectral distinguishability within pairs has a dra¬ 
matically deleterious effect on the attainable precision. 
Building on recent advances in spatially resolved single¬ 
photon detection [U [3Qti39] , we demonstrate that by con¬ 
trolling the input spatial structure of interfering photons 
and extracting complete spatial information at the detec¬ 
tion stage it is nevertheless possible to recover the sub- 
shot-noise precision. This confirms in proof-of-principle 
settings the feasibility of mode engineering techniques for 
quantum-enhanced interferometry. 


RESULTS 

Realistic two-photon interferometer. A generic 
two-photon Mach-Zehnder interferometer constructed 
with a pair of balanced 50/50 beam splitters and fed 
with photon pairs is shown schematically in Fig. Ba). 
The phase shift 0 between the interferometer arms modu¬ 
lates probabilities of detection events at the output ports 
which can be grouped into two types: either the pho¬ 
tons exit through different paths, producing a coinci¬ 
dence event between the detectors monitoring the ports, 
or both are found in the same output port leading to a 
double event. If the two photons are indistinguishable 
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FIG. 1. Two-photon interferometer, (a) In the multi- 
mode description, two paths are described by families of an¬ 
nihilation operators and subjected pairwise to a uni¬ 
tary map dependent on the phase shift 0. (b) Probabilities 
of coincidence pc and double count events pd at the inter¬ 
ferometer output are noticeably affected around ^ 7r/2 by 

imperfect indistinguishability (red dashed lines for the visi¬ 
bility V = 0.93) when compared with the ideal case (black 
solid lines). Grey dashed lines depict standard single-photon 
interference fringes, when only one input port is illuminated, 
(c) Residual distinguishability has a dramatic effect on the 
phase estimation uncertainty, shown with correspondingly 
coded lines, which diverges at 0 = 7r/2 for non-unit visibil¬ 
ity V. Engineering overlap in an additional degree of freedom 
of the interfering photons and implementing optimal measure¬ 
ment allows one to restore the sub-shot-noise precision over 
the entire operating range, shown with the purple solid line 
when the fraction of overlapping pairs is optimized individ¬ 
ually for each operating point. The blue solid line depicts 
an explicit strategy based on introducing a hxed transverse 
displacement between interferometer inputs and performing 
spatially resolved detection. These predictions are conhrmed 
by the experimentally determined estimation precision, de¬ 
picted as solid circles with two-sigma error bars. 


at the input, the first beam splitter generates a coher¬ 
ent superposition of both the photons in one or another 
arm of the interferometer, which is the simplest case of 
a NOON state providing sensitivity that approaches the 
Heisenberg limit pTHT^ . 

In the multimode description of the setup one intro¬ 
duces two sets of annihilation operators for upper path 
modes and lower path modes that are matched 
pairwise. Individual modes in each arm are mutually 
orthogonal, i.e. [d^,dj,] = The unitary 

map tJ{0) implemented by the interferometer between 
the input and the output ports transforms pairs of field 
operators labeled with the same index /i as 





( 1 ) 


Partial distinguishability of the interfering photons can 
be modeled by assuming that at the input the photon 
in the upper path occupies a certain mode di, while 
the lower-path photon is prepared in a combination of a 
matching mode bi and another orthogonal mode 62 with 
relative weights V and 1 — V, where V is the visibility pa¬ 
rameter specifying the fraction of indistinguishable pairs. 
To keep the notation concise we will write the complete 
two-photon state as pure 

IV^) = a\ {Vvb\ + Vl - Vbl) |vac), ( 2 ) 

and trace the final formulas over the index /i = 1,2 for 
the initial modes 6 ^. This is equivalent to taking from 
the start a reduced density matrix Va\b[ |vac) (vac|di 6 i + 
(1 —V)d| 62 |vac) (vac|di 62 . Here |vac) is the vacuum state 
of the entire multimode electromagnetic field satisfying 
d^jvac) = 6 ^ I vac) = 0 for any fi. Note that the above 
model includes the case of partly overlapping wavepack- 
ets constructed from a continuum of modes, wherein bi 
and 62 can be identified through the standard algebraic 
technique of Gram-Schmidt orthogonalization. 

The general transformation from Eq. 0 taken with 
/i = 1,2 implies the following expression for the proba¬ 
bility of a coincidence event 

PcW = 1 - ^(1 + V)sin2 0, (3) 


while the probability of a double event is Pd{0) = 1 — 
Pc{0). These formulas combine expressions for fully in¬ 
distinguishable and completely distinguishable pairs with 
respective probabilities V and 1 —V. The resulting fringes 
are depicted in Fig. 0 b) for V = 1 and 0.93. 

As a consequence of the Cramer-Rao bound m the 
minimum uncertainty of any unbiased phase estimate ob¬ 
tained from a measurement using N photon pairs around 
an operating point 0 is given by 


^pair _ 


(4) 
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where for standard photon counting at output ports 
of the interferometer the Fisher information is 

given by a sum of two terms contributed by coincidence 
and double events m, 


^pair(^) 


1 f dpcY 1 (dpdY 

Pe{d) \de J ^ pm \d0 ) 


(5) 


As a reference, we will take the uncertainty of ideal, shot- 
noise-limited phase measurement = l/v^2]V, when 

2N photons are sent individually to the interferometer. 
Our figure of merit will be the ratio e = of 

these two uncertainties, with e < 1 implying that sub¬ 
shot noise precision has been achieved. 

When the interferometer is fed with pairs of perfectly 
indistinguishable photons characterized by V = 1, we 
have 5 = I/V 2 independently of the operating point of 
the interferometer. It can be verified that around ^ = 0 
(equivalent to 0 = n) and 0 = nj2 the main contri¬ 
bution to Fisher information defined in Eq. comes 
respectively from double or coincidence events that oc¬ 
cur with vanishing probabilities when approaching these 
phase values. This is because in the ideal scenario even 
a small number of rare events provides a sound basis to 
infer the phase shift. Such a regime corresponds in stan¬ 
dard interferometry to dark-fringe operation used for ex¬ 
ample in gravitational wave detectors mm- 

As seen in Fig.[^c), the precision of phase estimation is 
affected dramatically by the non-ideal indistinguishabil- 
ity of photon pairs. In particular, statistical noise gener¬ 
ated by non-vanishing background of coincidence events 
effectively suppresses information about the phase shift 
that could be retrieved around 6> = 7r/2. We will refer to 
this operating point as the coincidence dark fringe. An 
analogous effect would be observed also at ^ = 0 if any 
mechanism generating spurious double events was incor¬ 
porated into calculations. 

Restoring quantum enhancement. The analysis 
presented above assumed that we have no access to the 
degree of freedom introducing partial distinguishability. 
For concreteness, we will consider this degree of freedom 
to be the spectral one, which means that all measure¬ 
ments performed on the photons are integrated in the 
frequency domain. Suppose now that we can fully con¬ 
trol and measure another, uncorrelated degree of freedom 
of the photon pairs sent to the interferometer. For the 
clarity of the argument, it will be convenient to use in 
this role the transverse spatial characteristics of the pho¬ 
tons. Let us consider a scenario when in addition to 
spectral distinguishability characterized by V we reduce 
the spatial overlap of the photons by preparing them in 
nonorthogonal spatial modes. As a result, even in the 
regime of perfect spectral indistinguishability only a frac¬ 
tion V of photon pairs would effectively overlap in space. 
To account for this scenario we will take the mode index 
to have two components /i = iy, where i = 1, 2 refers to 


the spectral degree of freedom, while y = R, L denotes 
two mutually orthogonal spatial modes. Using this nota¬ 
tion, the input state is described by an expression 

analogous to Eq. (§ with the following substitution of 
creation operators: 

bl ^ , i = 1,2 (6) 

and it explicitly reads 

\m = ab [vV(v^^b + VT^b\^) 

+ VTm2{VVblj^ + Vl^bl^)] |vac) (7) 

Note that both the spectral components bi and 62 have 
been subjected to the same spatial transformation. This 
is in accordance with our assumption that photon ma¬ 
nipulations cannot depend on the inaccessible spectral 
degree of freedom. 

Although the spectral and the spatial degrees of free¬ 
dom are treated on equal footing in Eq. Q, the crucial 
difference is the ability to measure the latter at the in¬ 
terferometer output. Therefore the phase shift 0 can be 
read out from the paths taken by the photons as well as 
their transverse spatial properties. To find the optimal 
strategy, we will resort to the concept of quantum Eisher 
information Fq(^) j42ll43] . which defines through an ex¬ 
pression analogous to Eq. Q the minimum uncertainty 
of a phase estimate inferred from the entire available 
characteristics of the physical system used for sensing. 
When the spatial degree of freedom for coincidence events 
is taken into account, the explicit expression for quan¬ 
tum Eisher information at the coincidence dark fringe 
0 = 7r/2, where the effects of spectral distinguishability 
are most severe, reads 

FQ{n/2) = 2^^. (8) 

Detailed derivation of this result is presented in Methods. 
Eor a given V, the maximum value of the above expres¬ 
sion is obtained for = (1 — \/l — V^)/V and reads 
Fg^^(7r/2) = 4(1 — Remarkably, this value 

gives sub-shot noise precision for any V > 0. In Eig.j^we 
compare the precision implied by numerically computed 
Fq{0) in the range 0 < ^ < tt with V optimized for an 
individual operating point to a scenario when no mode 
engineering has been attempted. A cross-section of these 
plots for V = 0.93 has been also shown in Eig. 00). It is 
seen that the singularity in precision around the coinci¬ 
dence dark fringe is removed and that the sub-shot noise 
operation is ensured across the entire range of 0. Note 
that this result is achieved without any postselection of 
two-photon detection events and no filtering or any other 
manipulation in the spectral domain has been applied. 
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FIG. 2. Enhancement of estimation precision. Relative 
enhancement £ = ^pair^^shot spatial overlap V = 

1 of the two input photons and (b) the overlap parameter V 
optimized individually for each given spectral visibility 0 < 
V < 1 and an operating point 0 < 0 < tt. The white area in 
(a) depicts the region e > 1 where sub-shot noise sensitivity 
is lost. It is seen that spatial mode engineering allows one 
to restore quantum enhancement across the entire parameter 
range. The uncertainty of the two-photon scheme is 

given by quantum Fisher information derived in Eq. (26). 


Measurement scheme. To elucidate the origin of 
the above effect it is instructive to analyze the oper¬ 
ation of the interferometer in the R/L basis of trans¬ 
verse spatial modes. Detecting both the photons in RR 
modes means that they overlapped spatially at the input 
and therefore underwent imperfect two-photon interfer¬ 
ence affected by non-unit spectral visibility V with fringes 
shown in Fig. Bb). On the other hand, combinations 
RL and LR at the output imply that the upper-path 
photon was initially in the mode R and the lower-path 
photon in the mode L. Consequently, both the photons 
propagated through the interferometer as independent 
particles exhibiting single-photon interference which for 
0 = TT 12 gives the steepest slope of interference fringes as 
seen in Fig. [^b). If the lower-path photon was prepared 


in a statistical mixture of R and L modes, only single¬ 
photon interference would provide information about the 
phase shift at 6> = 7r/2 and the shot-noise limit could not 
be surpassed. However, because the lower-path photon 
is sent into the interferometer in a superposition of the 
modes R and L, information from two- and one-photon 
interference can be combined in a coherent way through 
a suitable choice of the measurement basis at the interfer¬ 
ometer output. Strikingly, although neither one- nor two- 
photon interference used separately beats the shot noise 
limit itself, their coherent combination restores quantum 
enhancement of the measurement. 

As derived in Methods, the explicit form of the opti¬ 
mal measurement attaining quantum Fisher information 
at the coincidence dark fringe requires discrimination be¬ 
tween double events and two types of coincidence events 
corresponding to the following projections in the spatial 
degree of freedom: 


w = VttsI^^> + 5(\/FI±‘) i“> 

In Fig. we depict the resulting interference fringes for 
all three types of events when V is optimized for V = 93% 
at the operating point 0 = nj2. It is seen that coinci¬ 
dence events resolved in the ± basis indeed exhibit both 
one- and two-photon interference providing a rather steep 
slope at 6> = 7r/2, while their overall probability remains 
relatively low at this operating point. Combination of 
these features yields sub-shot noise sensitivity. 



FIG. 3. Optimized interference fringes. The probabili¬ 
ties of double pd (purple dashed line) and two types of coin¬ 
cidence pc± (purple solid lines) events corresponding to pro¬ 
jections in the optimal basis |±) defined in Eq. it with the 
overlap parameter R optimizing quantum Eisher information 
at ^ = 7r/2 for the visibility V = 0.93 in the inaccessible de¬ 
gree of freedom. Dashed red lines show standard two-photon 
interference fringes for the same visibility without engineering 
an additional degree of freedom. 
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In the limit V —)► 0 the optimal V approaches zero 
which means that at the input the two photons are nearly 
fully distinguishable in their spatial degree of freedom. In 
this regime the asymptotic expressions for the optimal 
measurement basis are |+) = \RL) and |—) = \LR) ^ i.e. 
we need to identify the origin of a photon that has ap¬ 
peared at a given output port. This is equivalent to real¬ 
izing single-photon interference twice, each time sending 
a photon into a different input port. The above anal¬ 
ysis explains why in the limit of zero visibility V = 0 
we recover the shot noise limit as evidenced by Fig. Hb). 
The described scheme is able to exploit non-classical two- 
photon interference for any V > 0 to achieve sub-shot 
noise operation, although unsurprisingly with a dimin¬ 
ishing quantum enhancement when V ^ 0. 

Spatially resolved detection. The measurement 
maximizing quantum Fisher information requires imple¬ 
mentation of rather exotic projections on two-photon su¬ 
perposition states given in Eq. This leads to the 

question about practical realization of the measurement 
achieving sub-shot noise precision at a given operating 
point. Fortunately, as we will now demonstrate, for 
generic spatial modes sub-shot noise precision can be re¬ 
stored just by measuring transverse positions of photons 
emerging from the interferometer. To discuss quanti¬ 
tatively this idea it will be convenient to resort to the 
paraxial approximation and to introduce two sets of op¬ 
erators ai{x) and bi{x) labeled with a continuous one 
dimensional transverse position x. As before, the index 
i = 1,2 labels spectral modes. Suppose now that the 
photons entering the interferometer along the upper and 
the lower path are prepared in spatial modes described 
by respective normalized profiles u{x) and v{x), 

aiR = J dxu{x)ai{x) 

VVblj^ + Vl - = J dxv{x)bi{x). (10) 

The overlap parameter now reads V = \ f dx {x)v{x)\‘^. 
The probability of detecting photons in two different out¬ 
put ports at positions x and x' is given by 

Pc{x,x'\e)= y] \{va.c\di{x)bj{x')U{d)\'tp'v)f (H) 

hi=i,2 


of photons detected in coincidence events. The probabil¬ 
ity distribution for this variable is given explicitly by 


1 cos" |e-(«-d)V4.^ 


Pc{^\d) = j dx'pcix'+ ^,x'\e) 

2ay/7r V 

+ sin" . (13) 

When the relative position of the two photons in coin¬ 
cidence events is available, Fisher information defined in 
Eq. ^ becomes enhanced by replacing the first term in 
the sum with the following integral over 


PciO) 


dpc 

de 


Fc{0) = I 


d^ 


dpc{^\0) \' 
Pc{^\d) \ 89 J 


(14) 

In Eig. I^c) we depict the estimation precision for the 
ratio dja = 1.64 which was used in the experiment de¬ 
scribed below. It is seen that the precision is brought to 
the sub-shot noise regime over the entire operating range. 

Experiment. To verify experimentally sub-shot noise 
phase sensitivity of the interferometric scheme described 
above we constructed an optical setup shown in Eig. 
The interferometer is fed with 800 nm photon pairs gener¬ 
ated via type-II spontaneous parametric down-conversion 
process, which are synchronized using a delayed line, spa¬ 
tially filtered through a single-mode fiber and delivered 
to the setup in two mutually orthogonal linear polariza¬ 
tions corresponding to the two input ports of the Mach- 
Zehnder interferometer. The photons emerging from the 
fiber were partly separated in space by inserting a 1.9 mm 
long calcite displacer D, which results in the displacement 
of d = 200 jam between the two ortogonally polarized 
output paths. The spatial modes can be modelled by 


gaussian functions defined in Eq. (12) with a = 122 /am. 


To ensure temporal stability, the interferometer trans¬ 
formation is implemented in the common-path configu¬ 
ration as a half-wave plate HWP, with the rotation angle 
equal to quadruple the phase shift 0 between the inter¬ 
ferometer arms, followed by a calcite beam separator S. 
The equivalence of this setup with the standard Mach- 
Zehnder interferometer is evidenced by decomposing the 
interferometer transformation introduced in Eq. § as 


where the summation over i = 1,2 stems from tracing 
over the spectral degree of freedom. 

As a concrete example we will take two gaussian modes 
of width a displaced by d, 

^a;) = vix)=uix-d) (12) 

VZTTCr^ 

that can be readily prepared by simple experimental 
means. In the phase estimation procedure we will exploit 
information contained in the relative position ^ = x — x' 



In the common-path configuration the diagonal matrix 
with entries 1 and —1 describes the half-wave plate in 
the coordinate system of its principal axes and the two 
outer matrices correspond to rotation to the laboratory 
reference frame by an angle d/4. 
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FIG. 4. Experimental setup. Interference between two orthogonally polarized photons a and h delivered by the single mode 
fiber SMF is realized in the common path configuration with the phase shift 0 corresponding to the quadrupled rotation angle of 
the half-wave plate HWP. The preceding calcite displacer D introduces transverse displacement between photon modes, which 
combined with spatially resolved detection restores the sub-shot-noise precision of phase estimation around the coincidence 
dark fringe. The output ports of the calcite beam separator S are mapped using a spherical lens SL and a cylindrical lens CL 
onto the intensified sCMOS camera detecting individual coincidence events with spatial resolution as shown in the inset. The 
upper right part of the figure depicts the spatial profiles of interferometer modes at the consecutive stages of the setup. 


We determined through an independent measurement 
(see Methods) the residual spectral distinguishability of 
the photons to yield V = 93%. Note that this value 
of visibility has been used in the example presented in 
Fig. [^b). The transverse distance between the two spa¬ 
tially separated modes corresponding to the output ports 
of the interferometer is 3.2 mm. The rear surface of the 
separator S is imaged by means of a spherical lens SL 
onto a single-photon-sensitive camera system (see Meth¬ 
ods) capable of spatially resolved detection of coincidence 
events [38l |39] . The profiles of spatial modes propagating 
in this case through consecutive stages of the setup are 
shown in the upper right part of Fig. 

We recorded spatially resolved two-photon detection 
events for three values of the phase shift around the dark 
fringe, 0 = 7r/2 — 0.1, 7r/2,7r/2 + 0.1, registering approx¬ 
imately 6 X 10^ events in each case. Experimentally 
observed spatial distributions Pc{x^x'\0) of coincidence 
events along with their marginals Pc{^\0) for ^ = x — x' 
are presented in Fig. It is seen that the joint posi¬ 
tion distributions are clearly sensitive to the phase shift 
in particular the sign of its deviation from 7r/2 can 
be unambiguously inferred from the asymmetry of the 
distribution with respect to the diagonal. 

To quantify information about the phase shift present 
in spatial distributions, we performed phase estimation 
from the actual experimental data and determined the 


estimation precision. In the preliminary step, we veri¬ 
fied the applicability of Eq. as a statistical model 
for the collected data assuming independently measured 
mode parameters a and d. We used the maximum- 
likelihood method to fit V and 0, obtaining for the 
three cases depicted in Eig.j^b) respective phase values 
0 = 1.47(2), 1.57(2), 1.69(2) and V = 0.93 which is in 
agreement with the visibility inferred from the indepen¬ 
dently measured Hong-Ou-Mandel dip. 

In order to determine the actual estimation precision, 
we divided data obtained for a given HWP setting into 
approx. 600 subsets of 10 two-photon detection events 
and estimated the value of the phase shift separately from 
each subset. The width of the resulting distribution of 
individual estimates can be used as a figure for the es¬ 
timation precision. The choice of the right estimation 
procedure needs some attention for small sizes of data 
sets. The maximum likelihood estimator is known to be 
asymptotically efficient [40] i.e. it saturates the Cramer- 
Rao bound in the asymptotic limit of infinitely many 
independent data samples. However, its application to 
small data sets is not justified owing to potential biased¬ 
ness. Therefore we used an estimator which is mani¬ 
festly unbiased for any data size and yields the precision 
given by Eq. at least in the vicinity of a given op¬ 
eration point. Specifically, for an experimentally mea¬ 
sured statistical frequency distribution /(^) of the rela- 
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FIG. 5. Spatial distributions of coincidence events, (a) Experimentally measured joint distributions pc{x^x\0) of 
detecting two photons at positions x and x' in two different output ports of the interferometer for photon pairs characterized 
by spectral indistinguishability V = 0.93 and three phase shifts 0 — 7r/2 — 12^7112 + 0.1. (b) Useful information about 

the phase shift can be extracted from the marginal histograms of the relative position between the photons ^ — x — x\ shown 
along with the fitted theoretical model Pc{^\0) used in the estimation procedure (solid lines). 


tive distance between the two photons, this estimator in 
the vicinity of = 7r/2 is explicitly given by [T7] 


^1/] = ^0 + 


=(0o) / 




fio dp,m 


Pc{C\do) de 


(16) 


e=eo 


where Fc{0) has been defined in Eq. (14). In practice, 
the integral in the above formula is discretized accord¬ 
ing to the width of the histogram bins. Applying this 
estimator to individual data subsets yielded a distribu¬ 
tion of phase estimates. The estimation precision was de¬ 
termined as the standard deviation of this distribution, 
shown in Fig. [^c) for the three phase shifts along with 
two-sigma error bars. The obtained values, clearly sit¬ 
uated below the shot-noise limit, demonstrate quantum 
enhanced operation of the interferometer despite partial 
spectral distinguishability of the input photons. 


DISCUSSION 

We analyzed operation of a two-photon Mach-Zehnder 
interferometer, which is one of the first [22] and 
most advanced technologically [6] examples of quantum- 
enhanced metrology, in a scenario which included two de¬ 
grees of freedom for interfering photons. It was assumed 
that one degree of freedom was inaccessible experimen¬ 
tally. Residual distinguishability of interfering photons 


in this degree of freedom has a dramatically deleterious 
effect on the precision of phase estimation around the 
operating point where coincidences at the interferome¬ 
ter output are suppressed owing to the Hong-Ou-Mandel 
effect, producing a coincidence dark fringe. We showed 
that exploiting another, completely uncorrelated degree 
of freedom over which one has full experimental control 
can mitigate this effect, restoring quantum-enhanced pre¬ 
cision in the entire operating range. This result is based 
on a rather subtle interplay between one- and two-photon 
interference. At the coincidence dark fringe imperfect 
two-photon interference alone does not provide any infor¬ 
mation about the phase shift, while one-photon interfer¬ 
ence obviously cannot surpass the shot-noise limit on its 
own. We demonstrated that combining coherently both 
types of interference through a suitably designed prepa¬ 
ration and measurement scheme in the second degree of 
freedom yields precision below the shot-noise limit. The 
feasibility of this approach was confirmed in an experi¬ 
ment using the transverse position of the photons as the 
controllable degree of freedom which could be measured 
with high resolution using a single photon sensitive cam¬ 
era system jSS] [39] . The presented strategy can be also 
applied to other scenarios involving two independent de¬ 
grees of freedom with different experimental accessibility, 
for example to mitigate effects of residual spatial distin¬ 
guishability when only area-integrating detection is avail- 
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FIG. 6. Three-photon scheme, (a) A Mach-Zehnder in¬ 
terferometer fed with a combination of two photons at one 
input port and a single photon at the other one. (b) Prob¬ 
abilities Pnn'i^) of detecting respectively n and n' photons 
at the interferometer outputs for perfectly indistinguishable 
photons V = 1 and the single photon exhibiting residual dis- 
tinguishability with respect to the other two mutually identi¬ 
cal photons, characterized by V = 0.93. (c) In the presence of 
residual distinguishability, the precision of phase estimation, 
shown with a dashed red line, deteriorates around operat¬ 
ing points where dominant contribution to Fisher information 
comes from the suppression of events nn — 21 or 12. This 
effect can be mitigated by introducing spatial displacement 
d = 1.45cr between interferometer inputs, resulting in quan¬ 
tum enhancement across the entire operating range (solid blue 
line). 


able, but photons can be suitably manipulated and mea¬ 
sured in the spectral domain. 

Beyond two-photon interferometry, one can consider a 
simple scenario involving a higher photon number when 
the input ports of a Mach-Zenhder interferometer are fed 
respectively with two photons and a single photon, as 


shown in Fig. |^a). If all three photons are indistin¬ 
guishable, counting photons at the interferometer out¬ 
puts yields fringes depicted in Fig. [^b) for which Fisher 
information = 7 irrespectively of the operating 

point. This is more than two-fold improvement over the 
shot noise limit for three photons used independently. 
When the single photon exhibits residual spectral distin¬ 
guishability with respect to the two photons (assumed 
to be identical) feeding the other port, minute changes 
of the shapes of the interference fringes have a dramatic 
effect around points ^0 = 2 arctan(l/\/2) and tt — Oq on 
a relative phase estimation uncertainty 5 = 
plotted in Fig.j^c), where A*^^^ is defined analogously to 
Eq. 0. According to Fig.j^b), this can be attributed to 
a non-vanishing background of events when two photons 
are detected at one output port and the remaining single 
photon at the other one. Preparing two photons in the 
same spatial mode u{x) and the single photon sent indi¬ 
vidually in a partly overlapping mode v{x) allows one to 
restore quantum enhancement as seen in Fig. [^c). In¬ 
terestingly, this effect occurs at different operating points 
compared to the two-photon case. 

A worthwhile candidate to analyze the benefits of mode 
engineering in a scalable multiphoton scenario may be 
the celebrated Holland-Burnett scheme [44] employing 
two Fock states with equal photon numbers. It is easy 
to verify that sub-shot noise precision at 0 = 7r/2 origi¬ 
nates from the suppression of odd photon number events 
at the interferometer outputs, which again is sensitive to 
residual distinguishability of input photons. A quanti¬ 
tative analysis of this scenario would require developing 
an efficient approach to deal with multimode multipho¬ 
ton states. On the other hand, it is known that in the 
presence of certain common imperfections, such as pho¬ 
ton loss, the ultimate precision follows asymptotically 
the shot-noise type scaling and the quantum enhance¬ 
ment has the form of a multiplicative factor which can 
be attained via a repetitive use of finite-size multiparticle 
superposition states m- Therefore results obtained for 
a fixed number of photons may also prove useful in the 
asymptotic limit for realistic scenarios. 

The presented results are based on the multimode 
description of a two-photon interferometer, which goes 
beyond the simplest models typically used to conceive 
quantum-enhanced measurement schemes. In practice, 
the applicability of quantum-enhanced techniques de¬ 
pends crucially on the ability to reduce decoherence ef¬ 
fects caused by noise and experimental imperfections. 
Results presented in this paper suggest that in addition 
to obvious attempts to suppress decoherence effects in 
interferometry by improving transmission of optical el¬ 
ements, stabilizing phase reference, etc., exploiting the 
multimode structure of quantum fields can help to to 
achieve the non-classical regime of operation. If the 
modal structure of the probes is carefully engineered at 
the input and suitably detected, such a strategy can of- 
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fer a noticeable improvement in precision even though 
decoherence effects due to another degree of freedom we 
do not have control over remain at the same level. An 
interesting question is whether analogous strategies can 
be generalized to quantum-enhanced interferometry with 
squeezed states of light [amiiis], boson sampling with 
linear multiport devices [THHH] , or perhaps benefit other 
quantum technologies such as optical quantum comput¬ 
ing [471148] . 

METHODS 

Output two-photon state. For two photons pre¬ 
pared initially in the state \ipT>) defined in Eq. 0 it will 
be convenient to write the output state as a sum of two 
components, 

u{e)\iJD) = \M0)) + \i’dm ( 17 ) 

Here |'0c(^)) is the conditional state describing coinci¬ 
dence events, when the two photons leave the interfer¬ 
ometer through different ports 

li’ciO)) = -[VWa\J)\j^cosd 

+ y/V{l-V) {alnblR cos^ f - alR^R f) 

+ cos^ I - sin^ |) 

+ y/{l-V){l-V){a\J)l^ cos^ |-4L^bsin^ f)]|vac), 

(18) 


while |'0d(6>)) corresponds to double events, when both 
the photons emerge at the same output port of the inter¬ 
ferometer: 

\M^)) = I sin 

+ \/(l — 'D)V{a\j^d\j^ — 

+ — 'D){1 — 

The overall probabilities for double and coincidence 
events are 

Pd{0) = {'ipd{0)\'ipd{0)) = i(l+ PV)sin^6> (20) 

and Pc{0) = 1 — Pd{0). Inserting V = 1 gives as a special 
case the result presented in Eq. 0. 

Let us note that all the terms in l^ipdi^)) exhibit iden¬ 
tical dependence on 0. Consequently, resolving double 
events with respect to the spatial degree of freedom can¬ 
not yield more information about the phase shift. There¬ 
fore we will focus our attention on coincidence events 
described by |'0c(^))- In order to analyze information 
about 0 when the spectral degree of freedom cannot be 
accessed, we will treat it formally as another subsystem 
ft, writing |vac) = |xx') (g) \ii')n, where i,i' = 1,2 

and x^x' = R^L. Tracing the two-photon state over 
the spectral subsystem yields the reduced density matrix 
Qc{^) = TrQ[|'0c(^)) ('0c(^)|] which written in the basis of 
spatial modes \RR), \RL) ,\LR), \LL) reads: 


^c(^) 


P[l-i(l + V)sin2 6>] 

yp(i — p)(cos^ i “ F ^) 
y/V{l-V){sm^ I - ir sin2 9) 
0 


yp(l - P)(cos4 I - iVsin^ 9) ^V{l-V){sm^ f - iVsin^ 9) 
(1 — P) cos^ I — j(l — P)Vsin^ 0 

-l{l-V)Vsm^9 (1-P)sin4| 

0 0 


0\ 

0 

0 

oj 

( 21 ) 


Fisher information. If the spatial degree of freedom 
of the photons at the interferometer output is projected 
onto the R/L basis, the complete statistics of measure¬ 
ment results is described by the diagonal elements of the 
density matrix gdO) for coincidence events and the col¬ 
lective probability Pd{0) for all double events. An easy 
calculation yields the corresponding Eisher information 


Fr/l{0) = 


2(1 + PV) - (P + 1)(V + 1) sin^ 9 
1 - i(V + l)sin2 6» 


( 22 ) 


At the dark fringe we have FjijL{7rl2) = 2{1 — V). This 
expression, which does not even reach the shot-noise 
level, can be understood intuitively: information about 
the phase shift is obtained only from mode-mismatched 
pairs, when the spatial mode or L at the output iden¬ 
tifies unambiguously the input port of a given photon. 


The factor 1 — P in F/^//^(7r/2) is the overall fraction 
of these events, while the constant 2 is contributed by 
single-photon interference exhibited by such pairs. 


Clearly, a measurement in the R/L basis neglects in¬ 
formation contained in the off-diagonal elements of the 
density matrix gdO). A measurement strategy that ex¬ 
ploits optimally gdO) is described by quantum Eisher 
information involving the symmetric logarithmic deriva¬ 
tive. Its calculation is simplified by switching to a new 
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basis for the spatial degree of freedom, 


\P) = ^{\RL)-\LR)), 

= + (23) 

The conditional density matrix gdO) takes the form 


{a\Qc{0)\l3) = {j3\Qc{9)\a) = i\/l -P2 cosO 

me{eW) = (i - (24) 

while all other elements involving jy) or (yj vanish. The 
symmetric logarithmic derivative Lc{0), given in general 
by the implicit formula 

^ = ^[L,{e)g,{e) + g,{e)U 0 )] ( 25 ) 

can now be easily found in the two-dimensional subspace 
spanned by ja) and \P). If double events are not re¬ 
solved in the spatial degree of freedom, quantum Fisher 
information can be written as 




FQ{e)=Tr{g,{e)[Ud)]^} + 


1 

Pd{0) 



^l-p2 + (l-|-PV)2cos2(9 
^1-PV + (1 + PV) cos2(9 


sin^ 0-\-2{l-\-W) cos^ 

(26) 


where Pd{0) is given by Eq. (20). Specializing the above 
expression to 9 = yields Eq. (|^, whereas its value 
for V optimized individually for a given spectral visibility 
0 < V < 1 and an operating point 0 < 6> < tt is presented 
in Fig. llJb). The symmetric logarithmic derivative has a 
simple off-diagonal form at ^ = 7r/2 


^<=(’^/2) = -2it:^(I«)(^I + l^)(«l)- (27) 

Quantum Eisher information is saturated by projecting 
the spatial degree of freedom onto the eigenstates of Lc{9) 
[121 [ 43 ] given explicitly for = 7r/2 by |±) = (|<a) ± 
|;d))/\/2. These states are nontrivial superpositions of 
photon pairs prepared in combinations of R^L spatial 
modes at two different output ports of the interferometer. 

Experimental details. In the experiment we used a 
photon pair source based on the Il-type SPEC process 
in a 5-mm long periodically poled KTP crystal (Raicol 
Crystals) pumped with 8 mW of 400 nm light from a 


continuous wave diode laser. The produced pairs are 
transmitted through a 3 nm EWHM interference filter, 
carefully synchronized in time using a delay line and spa¬ 
tially filtered using the single mode fiber. The gaussian- 
like spatial modes of the photons after the fiber have a 
fiat phase and the half-width a = 122 jam at 1/e height 
for the intensity distribution measured at the position of 
the camera system. The residual spectral distinguishabil- 
ity of the photons was determined from the depth of the 
Hong-Ou-Mandel dip scanned using the delay line and 
measured with standard avalanche photodiodes for the 
half-wave plate orientation corresponding to 6> = 7r/2. 


Our camera system [38l |39] begins with an image in- 
tensifier (Hamamatsu V7090D) where each detected pho¬ 
ton that induces a photoelectron emission produces a 
macroscopic charge avalanche resulting in a bright flash 
at the output phosphor screen. The flashes are subse¬ 
quently imaged with a relay lens onto a fast, low-noise 
6.5 pm X 6.5 pm pixel size sCMOS sensor (Andor Zyla) 
and recorded as approx. 25-px gaussian wide spots which 
can be easily discriminated from the low-noise back¬ 
ground. The central positions of the spots are retrieved 
from each captured frame with a subpixel resolution by 
a real-time software algorithm which provides full infor¬ 
mation about transverse coordinates of each registered 
coincidence event as illustrated in the inset of Eig. 

Eor the sake of simplicity we consider only the coordi¬ 
nate X in the horizontal plane of the setup and integrate 
the signals in the vertical direction. A cylindrical lens CL 
with / = 30 mm in front of the detector was used to re¬ 
duce the vertical size of the image, producing effectively 
a 700 px X 22 px stripe, which significantly decreases 
frame readout time and allows us to reach 7 kHz collec¬ 
tion rate of frames with exposure time 30 ns each. The 
overall quantum efficiency of the camera system is 23%. 


Three-photon scheme. A straightforward but te¬ 
dious calculation shows that for a three-photon input 
state of the form [VV[VRb\j^ + + 

vac) spatially resolved 
distributions for events 21 and 12 when the photons are 
split between the interferometer output ports into two 
and one are given by 


P2i{xi,X2]x'\0) = V\lsm0sm^{u{xi)v{x2) 

-^u{x2)v{xi))u{x') — COS^ ^u{xi)u{x2)v{x')\‘^ 
+(l-'^)[l sin^^sin^ |(|it(a:;i)|^|i;(a:;2)|^+|v(a:;i)|^|M(a:2)P) 
X \u{x')\'^ + cos^ ^\u{x{)\'^\u{x 2 )\'^\v{x')\'^^^ (28) 
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and 


Pi 2 {x; X2\0) = V |sin^ ^v{x)u{x[)u{x 2 ) 

— ^sinO cos ^u{x)(u{x[)v{x2) +u{x2)v{x[))\^ 

+ (1 - V)[isin2 6»cos2 ^\u{x)\'^{\u{x[)\'^\v{x'2)\'^ 

+ \u{x' 2 )\^\v{x[)\^) + sin® ^\v{x)\^\u{x[)\^\u{x' 2 )\^ 


(29) 


where we have made use of Eq. (10). Fisher information 


taking into account spatially resolved detection of events 
21 and 12 is given by 


F^^\0) = 


/ 

/ 


1 f dp30 \' 
P3o{0) \ dd ) 
dxi dx 2 dx' 

P2l{xi,X2]x'\0) 

dxdx[dx2 

Pi2{x;x[,x^2\^) 


+ 


1 


f dpQ3 \' 
PO3{0) \ d0 ) 


■^P2i{xi,X2;x'\0) 

d 

■^Pi2ix;x[,X2\0) 


(30) 


Assuming Gaussian spatial modes introduced in Eq. (12), 


the above expression was optimized over the displace¬ 
ment d for Oq = 2 arctan(l/v/2) and V = 93%. The 
relative uncertainty of a phase estimate for the obtained 
value d = 1.45cr has been depicted for the entire range 
0 < < TT in Fig. [^c). 
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